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Open Fabry-Perot microcavities represent a promising route for achieving a quantum electrody¬ 
namics (cavity-QED) platform with diamond-based emitters. In particular, they offer the oppor¬ 
tunity to introduce high purity, minimally fabricated material into a tunable, high quality factor 
optical resonator. Here, we demonstrate a fiber-based microcavity incorporating a thick (> 10 
pm) diamond membrane with a finesse of 17,000, corresponding to a quality factor Q ~ 10®. Such 
minimally fabricated, thick samples can contain optically stable emitters similar to those found in 
bulk diamond. We observe modified microcavity spectra in the presence of the membrane, and 
develop analytic and numerical models to describe the effect of the membrane on cavity modes, 
including loss and coupling to higher-order transverse modes. We estimate that a Purcell enhance¬ 
ment of approximately 20 should be possible for emitters within the diamond in this device, and 
provide evidence that better diamond surface treatments and mirror coatings could increase this 
value to 200 in a realistic system. 

PACS numbers: 42.81.Wg, 42.25.Hz, 42.60.Da, 42.50.-p 


I. INTRODUCTION 

The field of diamond photonics has seen tremendous 
growth over the last decade mm, spurred by new ap¬ 
plications of optically active defect centers in metrol- 
ogy 121 m and quantum information science [5j [6]. In 
particular, the nitrogen-vacancy (NV) defect [7] exhibits 
long spin coherence times and narrow optical transitions 
favorable for realizing a solid-state cavity-QED system. 
In pursuit of this goal, much progress has been made 
in fabricating low-mode-volume cavities in diamond it¬ 
self O ISUTO]. and Purcell enhancement of the NV zero 
phonon line as large as 70 has been observed in a diamond 
photonic crystal cavity m- A complementary strategy 
is to confine the defect in an open Fabry-Perot micro¬ 
cavity [nun], which provides in-situ tunability and the 
possibility for very narrow cavity linewidths. Recently, 
three groups have observed coupling between an open 
cavity and an NV center in a nanocrystal [HHii]. 

A central challenge for diamond photonics is desta¬ 
bilization of defect optical transitions in close proxim¬ 
ity to surfaces, especially for defects in nanocrystals or 
in nanofabricated devices m- For example, the afore¬ 
mentioned diamond photonic crystal cavity achieved its 
high Purcell factor at the expense of spectral diffusion of 
many GHz m, far in excess of the near lifetime-limited 
linewidths of 13 MHz that can be observed in some type 
Ila samples m- While recent advances in surface treat¬ 
ments [20] and fabrication [21] hold promise for realiz¬ 
ing optimal NV properties in nanophotonic structures, 
narrow linewidths are most reliably obtained microns 
into bulk diamond. The open cavities discussed here 
can accommodate the larger mode volumes required for 
such microns-thick samples: both their mode volume V 
and quality factor Q increase approximately linearly with 
length, so that the Purcell enhancement Fp depends only 
on the mirror finesse F and the ratio of the cavity waist 
wq to the resonant wavelength A: Fp oc FX^jw^ [12]. 


In addition, their linewidths are comparatively narrow 
and can be tuned over a wide range in-situ via the 
cavity length, potentially allowing exploration of spin- 
dependent coupling between an NV center and the cav¬ 
ity. Finally, by adjusting the positions of the mirrors, 
one can optimize the cavity mode spatial overlap with 
the emitter. 

In principle, it is straightforward to incorporate a 
microns-thick membrane into an open cavity. However, 
for high finesse F > 10^ cavities, losses at the 100 ppm 
level are important. Absorption and scattering must be 
minimized, and changes in the cavity mode induced by 
the dielectric interface must be considered. Here, we 
demonstrate that a fiber-based microcavity can main¬ 
tain high finesse F ^ 17,000 (quality factor Q ^ 10^) 
while incorporating a > 10 pm thick diamond membrane 
compatible with high stability defect centers. We fur¬ 
ther develop a theoretical description for the longitudinal 
modes (including diffraction effects), and perturbatively 
estimate the mixing between transverse modes induced 
by the membrane. Based on our measurements, we pre¬ 
dict our device should be capable of enhancing the NV 
zero phonon line by a factor of approximately 20. 


II. THE FIBER CAVITY DEVICE 

We work with a fiber-based Fabry-Perot microcav¬ 
ity [HI [13] in a geometry similar to those used to study 
quantum dots [22] and molecules [23]. The microcavity 
system (Fig. [^) consists of a concave mirror on the tip 
of a single mode optical fiber, and a macroscopic fiat mir¬ 
ror to which we bond the diamond membrane. Compared 
with traditional optics, these fiber-based cavities offer ad¬ 
vantages in stable alignment and efficient coupling to the 
single mode propagating in the fiber [12] . 

The fiber mirror substrate is fabricated using a CO 2 
laser ablation process [12]. The ablation laser induces 
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FIG. 1. (Color online) (a) A schematic of the microcavity sys¬ 
tem illustrating the cavity length (L), mirror radius of curva¬ 
ture (R), diamond membrane thickness (t^), and coupling ef¬ 
ficiency (e) between the fiber and cavity modes, (b) An inter¬ 
ferometric image of an ablated fiber, where each subsequent 
dark fringe corresponds to a change in depth of ~ 253 nm. (c) 
A microscope camera image of the fiber cavity mirror and re¬ 
flection seen against the diamond membrane. Dark regions in 
the upper corners are due to partial etching of the mirror be¬ 
neath the membrane, and some contamination of the diamond 
surface is also visible, (d) Measured transmission (solid line) 
and reflection (dashed line) curves for the bare cavity, and 
membrane-in-cavity configurations, normalized to the peak 
reflected power. These particular measurements correspond 
to finesses of JT ~ 37,000 ± 1000 and JT ~ 15,900 ± 900 for 
the bare cavity and membrane-in-cavity respectively. 


evaporation and melting of the glass on the fiber tip 1211 , 
resulting in an approximately Gaussian-shaped dimple 
(Fig. [^) with an extremely low surface roughness of 
< 0.2 nm-rms, as measured with an atomic force mi¬ 
croscope. By imaging the fiber core during ablation 
alignment, we achieve a placement repeatability of 0.5 
|LLm for the ablation spot. The fiber used in our exper¬ 
iments is measured to have a power coupling efficiency 
to the cavity mode of = 48 ± 4% (limited by mirror 
absorption as well as cavity and ablation misalignment, 
see Appendix |^, and an effective radius of curvature 

= 61.0 ± 1.4 |LLm (see Appendix [B|) . 

The second mirror substrate is a superpolished macro¬ 
scopic mirror flat with surface roughness below 0.1 nm- 
rms. The flat and fiber mirrors are coated with a di¬ 
electric mirror stack (LASEROPTIK) specified to have 
a transmission of 70 ± 10 ppm and < 24 ppm losses at 
A = 637 nm; the theoretical finesse of the stack design 
is 53,100. The fiber (stripped of all polymer jacketing) 
and flat mirrors were both annealed at 300° C for five 
hours under atmospheric conditions to reduce losses in 
the coatings m- 

The diamond membrane is fabricated from a (100)- 
cut electronic grade single crystal diamond plate. The 


bulk diamond was laser cut laterally, producing 20 ± 10 
|Lim thick diamond membranes polished to a surface 
roughness of approximately 5 nm-rms. One of the re¬ 
sulting membranes was cleaned in a piranha solution 
and bonded with Van der Waals forces to a silicon car¬ 
rier wafer, and approximately 2 \xm was etched from 
the membrane surface using an ArCl 2 inductively cou¬ 
pled plasma reactive ion etching (ICP RIE) recipe [26l - 
1^ . The etching process reduced the surface roughness 
to < 0.2 nm-rms (measured over one optical wavelength 
squared). The membrane was then removed from the car¬ 
rier, similarly cleaned in piranha, etched on the other side 
(again removing 2 fim with ArCl 2 ), and finally bonded to 
the macroscopic mirror flat. Einally, a third ArCl 2 etch 
was performed to thin the membrane to approximately 
10 |LLm. 

To assemble the cavity, the mirror flat is fixed to a tip- 
tilt mount, while the fiber is clamped to a 3-axis man¬ 
ual and piezo stage. The tip-tilt mount enables angular 
alignment of the cavity mode, while the 3-axis stage al¬ 
lows for study of different regions of the membrane as 
well as precision control of the cavity length. Eigure 
shows an image of the assembled device. 


III. EFFECT OF THE MEMBRANE ON 

CAVITY MODES 

A. Cavity finesse 

Introducing the membrane into the cavity affects the 
linewidth of its resonances. To determine the cavity fi¬ 
nesse, we scan the position of the fiber mirror while mon¬ 
itoring the cavity’s transmission and reflection at a fixed 
wavelength near 637 nm (provided by a tunable diode 
laser). In this case, we define the finesse as the ratio of 
the free spectral range (ESR « A/2) to the resonance full 
width half maximum (EWHM) measured as a function 
of the length of the cavity, and denote it by JV- 

_ ESR in length 

^ “ EWHM in length ^ ^ 

Note that for the membrane-in-cavity system, this is not 
necessarily the same as finesse obtained by measuring the 
resonance spacing and linewidth as a function of laser 
frequency. 

We measure the finesse by first performing a long scan 
of the cavity length to observe the resonance spacing as 
a function of the voltage applied to the piezo stage. By 
scanning the length over about 20 \im (roughly 60 ESR), 
we can fit the observed resonances to extract the free 
spectral range and calibrate the piezo stage nonlinearity. 
Subsequent voltage scans over shorter length ranges (0.6 
|LLm) provide high resolution data for extracting the cav¬ 
ity linewidth. At each position of interest, we measure 
64 transmission and reflection peak data sets to gather 
statistics on the cavity linewidth. This procedure is fol¬ 
lowed for all measurements of JV presented in Eigs. 1-4. 
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We characterize the transmission and reflection curves 
for the TEMoo fundamental mode using the same 
fiber mirror for an empty or “bare” cavity, and for a 
membrane-in-cavity, as illustrated by sample data sets in 
Fig.[T]i. These measurements are corrected for calibrated 
losses in the measurement apparatus (outside of the cav¬ 
ity), and are normalized to the peak reflected power. To 
determine cavity linewidth, we fit the transmission and 
reflection data sets to Lorentzian and Fano lineshapes re¬ 
spectively; in practice, the reflection signal gives better 
signal to noise and was used to calculate cavity finesse. 
The asymmetric resonances seen in reflection can arise 
from a slight displacement of the ablation dimple from 
the fiber core (see Appendix 0 . 

As discussed below, the observed finesse varies with the 
length of the cavity, the transverse location on the mem¬ 
brane, and the frequency of the laser. At best, we observe 
a peak finesse of ~ 37, 000 for the bare cavity and 
Tl ~ 17, 000 for the membrane-in-cavity setup. When 
we observe different locations on the flat mirror’s surface, 
we find the finesse for the empty cavity typically fluctu¬ 
ates by a few thousand, most likely due to surface con¬ 
tamination or spatially-varying surface roughness. With 
the diamond present, however, the finesse fluctuates by a 
much larger factor, with no observable cavity resonances 
in many locations. At first glance one might presume 
these fluctuations arise from similar physics, i.e. rough¬ 
ness, contamination, or even crystal defects in the di¬ 
amond itself. However, as discussed below, such large 
finesse fluctuations are primarily caused by spatial vari¬ 
ations in the diamond layer thickness, which affects the 
cavity mode structure in an important and predictable 
way. 


B. Mode structure 

We characterize the cavity mode structure by illumi¬ 
nating the flat mirror with a broadband LED source, and 
measuring the spectrum of the light transmitted into the 
fiber with a grating spectrometer. By gathering data as 
a function of cavity length, we observe the evolution of 
multiple longitudinal and transverse modes (see Fig.[^). 

The measured white light spectrum exhibits a canted 
periodic structure that is markedly different from the be¬ 
havior of a bare cavity. These features can be quan¬ 
titatively reproduced by a simple one-dimensional (ID) 
model. We consider lossless mirrors at each end of the 
cavity, with a 180 degree phase shift on reflection (facing 
the cavity) to approximate the dielectric mirror stack ter¬ 
minated at the high index material (Ta 205 in this case). 
Between the mirrors are a slab of diamond of thickness 
td and index Ud^ and a layer of air with thickness L — td 
and index riair = 1 (see Fig. [^). In the limit of perfect 
mirrors, the resonant frequencies u are given by solutions 


to the transcendental equation 

(1 + Ud) sin {L + tdirid - 1))) 

= {1 - Ud) sin {L - tdiud + 1))). (2) 

Note that while the resonances occur regularly every 
c/2z/ as the length of the cavity shifts, the variation 
with frequency is less straightforward. For long cavi¬ 
ties, encompassing many nodes, Eq. can be approx¬ 
imated by writing u in terms of its deviation from an 
integer multiple m of the average free spectral range, 
u = Su ^ mcl{2(L + {ud — l)^d)), and neglecting 6u in 
the RHS of Eq. [29], yielding: 

- (-1)" (gjl .in [ W 

Fitting Eq. [^to the fundamental mode frequencies in 
the cavity spectrum results in an estimated membrane 
thickness of td = 10.b ±0.2 jam and the cavity lengths 
given in the x-axis of Fig. [^, where the fit results are 
shown by the dashed lines. Note that we fit resonances 
over the full 20 jam range of the stage (not shown) to pro¬ 
duce these estimates, and included a cubic nonlinearity 
in the piezo stage response; the region displayed in the 
figure is representative of the goodness of fit. The fit de¬ 
viations arise because the model neglects the transverse 
Gaussian field profile and Guoy phase of the cavity mode, 
which can lead to errors in the estimated length of up to 
half a FSR (A/4) [30]. The fit also allows us to determine 
the cavity length during transmission and reflection mea¬ 
surements, albeit with an increased uncertainty (roughly 
±0.3 fim) [3l]. 

The model also provides some intuition about the sys¬ 
tem. If the membrane-air interface were perfectly reflec¬ 
tive (i.e. Ud oo), it would divide the cavity into two, 
and the normal modes would separate into “diamond” 
modes and “air” modes, wherein the field is entirely lo¬ 
calized in either the diamond or air, respectively. Since 
the diamond thickness is fixed, the diamond mode fre¬ 
quencies (horizontal lines in Fig. 2b) would not depend 
on the longitudinal position of the fiber mirror, while the 
air modes would decrease in frequency as the air gap in¬ 
creases in length (slanted lines in Fig. 2b). Indeed the 
frequency spacings for these modes would reveal the di¬ 
amond thickness and cavity length: Audi = 2 ^ 774 ’ 
Ai^air = 2 {L-ta) diamond and air modes respec¬ 

tively. With finite Ud^ these modes are coupled to one 
another, leading to the large avoided crossings observed 
in the spectrum; “diamond-like” modes have a shallow 
slope, while the “air-like” modes have a steeper slope. 
This behavior is very similar to that of a membrane-in- 
the-middle system [29] , where the air-diamond dielectric 
interface plays the role of a weakly reflective, vanishingly 
thin membrane. 
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FIG. 2. (Color online) (a) Cavity spectrum obtained by coupling a broadband LED through the flat mirror, and scanning the 
length of the cavity. An approximate analytic ht (Eq. ^ to the resonances is overlaid. The error in the x-axis calibration is 
±0.2/im. (b) A numerical simulation of the normalized hnesse Tl of the cavity resonances including the Guoy phase and mirror 
stack. Lines indicating the cavity resonances associated with air and diamond regions are overlaid to illustrate the avoided 
crossings, (c) Einesse vs. laser wavelength measured for three different regions of the diamond membrane, corresponding to 
different membrane thicknesses. The raw hnesse data was binned using the Ereedman-Diaconis rule to show the underlying 
distribution, and the binned data is shown through opacity. The mean hnesse values are plotted as lines, (d) Simulated hnesse 
vs. laser wavelength for a lossless cavity. The absence of scattering losses leads to high hnesse over a much larger wavelength 
range. 


Our analysis above focused on the fundamental mode, 
and indeed, for an ideal spherical ablation dimple, light 
from the hber core should couple primarily to the Gaus¬ 
sian TEMoo mode. Nevertheless, some higher order 
modes are also visible in the spectrum. Similar features 
observed in the white-light spectrum for the bare cavity 
are used to extract the effective radius of curvature of the 
hber mirror [R = 61.0 ± 1.4 qm), as noted in section [n| 
and detailed in Appendix [B| 

In addition to the analytic ID model used to hnd the 
cavity resonance frequencies, we developed a numerical 
three-dimensional model for the cavity modes that incor¬ 
porates wavefront curvature within the cavity and the full 
dielectric mirror stack. This model includes an approxi¬ 
mation that the air-diamond interface follows the curva¬ 
ture of the Gaussian wavefronts, in order to prevent cou¬ 
pling between transverse modes via refraction. Figure 
shows simulated cavity resonances in the absence of any 
loss, calculated over the same length and frequency range 
as the white light transmission data. The cavity param¬ 
eters used are those extracted from the ht in Fig.|^, and 
full calculation details are provided in Appendix]^ The 
color of each data point shows the calculated value of 
Rl normalized to the naive hnesse estimate of where 
T is transmission per pass of one mirror, and all other 
loss processes are neglected. The highest hnesse values 
are obtained when the laser frequency is tuned to an air¬ 


like mode, approaching the naive estimate. Conversely, 
if the laser frequency is tuned to a diamond-like mode, 
the measured hnesse will be consistently lower than ex¬ 
pected. 

We emphasize that our models thus far assume no 
losses, meaning the aforementioned hnesse huctuations 
arise entirely from interference effects. The hnesse lim¬ 
itations can be understood by considering the effect of 
attaching a diamond membrane to the hat mirror. The 
dielectric coatings used for our mirrors are terminated 
with a high index material, and are optimized for use in 
air. Diamond has a high index of refraction {ud = 2.417), 
which effectively lowers the rehectivity of the hat mir¬ 
ror, corresponding to a decrease in hnesse for modes in 
which the electric held is more conhned to the membrane. 
Quantitatively, the lossless ID analytic model predicts 
that the hnesse of the diamond-like modes is reduced by 
a factor of 2/(n^ + 1) « 0.3. Conversely, if low-index- 
terminated mirrors were used, the diamond-like modes 
would exhibit the naive hnesse while the air-like modes 
would have hnesse reduced by 2/(l/n^ + l) ^0.6. 

A central prediction of this calculation is that the mode 
structure can cause dramatic hnesse variations with laser 
wavelength. Moreover, because the mode structure shifts 
with td^ peak hnesse values should occur at different fre¬ 
quencies for different membrane thicknesses. Figure 
shows the measured hnesse as a function of laser 
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wavelength for three regions on the diamond membrane 
with marginally different thicknesses [32]. The square 
data points correspond to raw finesse data (binned us¬ 
ing the Freedman-Diaconis rule to show the underlying 
distribution) where opacity illustrates measurement fre¬ 
quency. The mean finesse is plotted with a line. For each 
region, the peak finesse occurs at some wavelength corre¬ 
sponding to an air-like mode. The finesse decreases as the 
laser is tuned away from this wavelength, and the electric 
field becomes more localized in the diamond membrane. 

While the qualitative features of our data in Fig. [^ 
are similar to the lossless model predictions, the drop in 
finesse is notably larger and steeper. For comparison. 
Fig. Ill shows the numerically simulated finesse for the 
ideal lossless system, which exhibits much more gradual 
variations. As discussed quantitatively below, the dis¬ 
crepancy can only be explained by including loss primar¬ 
ily at the air-diamond interface, such as scattering from 
roughness or contamination. The sharp wavelength de¬ 
pendence again arises from interferometric effects: when 
there is an electric field node at the air-diamond inter¬ 
face, field driven surface losses are strongly suppressed. 
In this geometry, a node appears at the air-diamond sur¬ 
face only for the air-like mode, providing a mechanism 
for the sharp finesse peaks in Fig. H- 

To quantitatively understand the effects of loss, we 
add different absorption and scattering mechanisms to 
the numerical transfer matrix model described in Ap¬ 
pendix 0 We consider loss in the mirrors, loss caused 
by scattering at the diamond interfaces, and absorption 
in the diamond. Loss inside the mirrors and diamond 
is modeled by adding complex components to the refrac¬ 
tive indices of the layers. Scattering by surface roughness 
of the diamond membrane is added by adjusting the in¬ 
terface reflection and transmission coefficients according 
to [331134] 


_ (0) -2(27rcrni/A)^ 

f.. _ ^(0) -(l/2)(27rcr(ni-nj)/A)2 


( 4 ) 

( 5 ) 


where rij (tij) is the amplitude reflection (transmission) 
coefficient going from material of index rii into mate¬ 
rial of index n^, a is the rms surface roughness, A is the 

wavelength in vacuum, and rand are the lossless 
Fresnel coefficients. These reflection and transmission co¬ 
efficients are used in the transfer matrix describing each 
diamond surface. To quantitatively compare the effects 
of each individual source of loss, we increase its strength 
sufficiently to bring the peak finesse J'l down to our ob¬ 
served value of 17,000 (while holding other sources of 
loss at zero). We then calculate the cavity modes and 
linewidths for the cavity parameters extracted from the 
fit in Fig. ^ over the cavity lengths illustrated in Fig.[^. 

Figure ^ shows Jb: as a function of wavelength, as 
predicted by several loss models. As noted earlier, scat¬ 
tering at the air-diamond interface behaves qualitatively 
differently from the other loss models, and most closely 
approaches the features we observe in Fig. H- Because 


there is always a node at the high-index-terminated mir¬ 
ror surface, scattering from the diamond-mirror interface 
does not produce such sharp features. 

Figuresalso includes simulations using our best esti¬ 
mate for the specific losses in our system. Enough mirror 
loss was added to bring the finesse down to 37,000 (peak 
finesse measured for the bare cavity), and we set the 
diamond-mirror interface roughness to 0.19 nm-rms (as 
measured for similar samples). To match the features in 
Fig.[^, we added sufficient scattering at the air-diamond 
interface to produce a peak in Tc (see Fig. |^) with a 
FWHM of 1.14 nm (the linewidth of the central peak in 
Fig. 1^). Finally, absorptive loss was added to bring the 
peak Tl value down to 17,000. Notably, a very large 
air-diamond surface roughness (a = 3.5 nm-rms) was re¬ 
quired to reproduce the features of Fig. H- This rough¬ 
ness is far larger than values <0.2 nm-rms measured on 
diamond samples etched by ArCl 2 , and indicates some 
additional surface scattering or contamination is likely 
to blame. 

After optical measurements were concluded, AFM 
measurements of the membrane revealed surface rough¬ 
ness of ^ 1 nm-rms in the regions of interest (measured 
over one optical wavelength squared). The increased 
roughness was likely caused by the third ArCl 2 etch 
(while bonded to the mirror), which produced notice¬ 
able surface damage in some areas of the membrane (see 
Fig. [^); the regions used in these experiments appeared 
unaffected, but in fact suffered roughening. Surface ab¬ 
sorption caused by contamination could also be present, 
and would behave similarly in our model, exhibiting the 
same sensitivity to a field node at the air-diamond inter¬ 
face. It is therefore likely that reduced losses could be ob¬ 
tained with better surface preparation. Based on our sim¬ 
ulations, state-of-the-art surface roughness a = 0.19 nm- 
rms is compatible with finesse > 10^ for air-like modes 
and > 50, 000 for diamond-like modes. 


IV. ESTIMATION OF PURCELL 
ENHANCEMENT 

A figure of merit for cavity systems is the Purcell fac¬ 
tor Fp, which describes the spontaneous emission en¬ 
hancement into the cavity mode for an optimally placed 
single emitter. For a cavity with varying refractive in¬ 
dex [351 ESI, 

p _ l-^maxP _ 

^ 47r^ rid fff (i*) (i*) * 

Here, Aiy is the cavity linewidth in frequency, A is the 
resonant wavelength, E(r) and n(r) are the electric field 
and index of refraction within the cavity, and F^ax is 
the electric field at the emitter in the diamond, assuming 
perfect emitter orientation and location. Notably, the 
Purcell factor depends on the linewidth in frequency, not 
length. As we are not currently able to directly measure 
the cavity spatial mode and Az/, we use our model and 
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FIG. 3. (Color online) Numerical cavity simulations in¬ 
cluding wavefront curvature and loss as a function of reso¬ 
nance wavelength, (a) Simulated cavity transmission, illus¬ 
trating the cavity lengths for which the simulation was per¬ 
formed (dashed line). Steeper-sloped regions correspond to 
diamond-like modes, while shallower-slopes correspond to air¬ 
like modes.(b) Simulated finesse J^l for different sources of 
loss in the membrane-in-cavity system. For each simulation, 
enough loss was added to bring the measured finesse in length 
down to 17,000. (c) Simulated Purcell enhancement for the 
same models as in (b). Note that the legend at bottom applies 
to both (b) and (c). 


measurements of J'l to provide a theoretical estimate of 


the Purcell enhancement available in this cavity geome¬ 
try. 

In our analytic ID model, it is straightforward to cal¬ 
culate the mode integrals of Eq. and the linewidth Au 
in terms of the mirror finesse (a function of only the 
lossless mirror reflectivity). This yields a simple result in 
the limit of large 


F^^'> = T 


6A^ 


Fy) = T 


12A2 

{nl + ridWw^ ’ 


(7) 

(8) 


where Fp'^'^ (Fp^^) is the Purcell enhancement for air¬ 
like (diamond-like) modes and wq is the 1/e^ intensity 
radius of the cavity waist. Note that in the limit ^ 1, 
these match what one would obtain from the standard 
Purcell formula (3A^/47r^)((5/E) with a mode volume of 
V = i7r/4)wlL and Q = 2LJ’/A. 

The two types of modes have different Purcell factors 
because they have different vacuum electric field max¬ 
ima in the diamond and different cavity linewidths in 
frequency. The variation in Az^ has contributions from 
the reduced reflectivity of the fiat mirror (due to the di¬ 
amond layer) as well as the relative round-trip times of 
the diamond and air half-cavities. Such effects are simi¬ 
lar to finesse oscillations observed in optomechanical sys¬ 
tems [37]. Remarkably, in the high finesse limit of the 
lossless ID model, the length dependence of the vacuum 
electric field maximum in diamond precisely cancels the 
length dependence of Az/, yielding the simple expressions 
above. 

For lossless systems, matches the peak value of 
and one might be tempted to use Eqs. [7]|^ with our ob¬ 
served peak finesse and cavity geometry to determine the 
Purcell enhancement of our device. Such a calculation 
(using T = 17,000, reo=2.2 qm, A = 637 nm) would pre¬ 
dict ^ 20 and ^ 33. However, adding in loss 

does not simply reduce F: the location of the loss (in 
diamond or air) will affect the modes differently, and in 
general we find that using Eqs. [7][^ with F = Max(J^i,) 
overestimates the best Purcell enhancement for realistic 
systems where loss is associated with the diamond. 

Eigurej^ shows the Purcell factor calculated using the 
numerical model with the different loss mechanisms de¬ 
scribed in the previous section. While mirror absorption 
produces similar results to the predictions of Eqs. [7]|^ 
qualitatively distinct behavior appears from the surface 
losses that likely limit our system. In particular, we pre¬ 
dict a maximum Purcell enhancement of approximately 
20 for our current device geometry. However, our anal¬ 
ysis also suggests that significant improvements can be 
obtained. Eor example, if surface losses can be limited to 
the observed roughness after ArCl 2 etching (< 0.2 nm- 
rms), and higher reflectivity mirror coatings are used, a 
cavity finesse of 50, 000 can be maintained even with an 
antinode at the air-diamond interface. Using a 30 jam 
radius-of-curvature mirror (attainable in our laser abla- 
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tion setup), a (lll)-oriented 5 /iui thick membrane, and 
a cavity length of 10 /im, a maximum Purcell factor of 
around 200 could be reached. Such a cavity would also 
couple efficiently to the fiber mode (> 90% with perfect 
alignment and low-loss mirrors [38]) and have a linewidth 
^300 MHz, which is large enough to accommodate mi¬ 
nor spectral diffusion but small enough to resolve the 
excited-state structure of the NV center. 


V. FINESSE CHANGES WITH CAVITY 
LENGTH 

Beyond the absorption and scattering processes con¬ 
sidered above, a thick diamond membrane could also in¬ 
duce an additional, potentially important source of loss: 
mixing between transverse modes of the cavity. Our nu¬ 
merical model has assumed that the air-diamond inter¬ 
face follows the spherical wavefront of the cavity mode, 
allowing description of the cavity eigenstates in terms 
of two Gaussian beams in the diamond and air regions. 
The real planar interface, however, deviates from this re¬ 
quirement, and can thereby couple the TEMqo mode into 
higher-order Hermite-Gaussian modes. Because high- 
order transverse modes have a larger spatial extent, this 
mechanism could induce additional losses caused by clip¬ 
ping at the small fiber mirror. This type of loss would be¬ 
have differently than those discussed previously because 
it would depend on the length of the cavity, with greater 
losses expected when higher-order modes approach de¬ 
generacy with the fundamental. 

To estimate such losses, we apply non-degenerate per¬ 
turbation theory (see Appendix 0 to calculate the 
first-order eigenstates of the membrane-in-cavity. The 
fraction of those eigenstates clipped at the fiber mir¬ 
ror can then be calculated to determine the loss per 
round trip. We begin with zero-order modes from the 
three-dimensional model discussed in Appendix We 
then treat the volume of diamond between the curved, 
wavefront-matching surface and the true flat interface as 
the perturbative volume. The first order correction to 
the eigenstate is given by [39] : 

I^m 1^00 

where the eigenstate is ^ 0oo + 4^00 is a TEMqo 

zero-order cavity mode derived from our model, 'iprn is 
the zero-order mode (including all longitudinal and 
transverse mode numbers). Km = (^) contains the 
eigenfrequency uUm for the mode, and the integral is 
taken over the volume of the perturbation. The overlap 
integral couples only even order transverse modes, and 
falls off quickly with transverse mode number. In prac¬ 
tice, we have included transverse mode numbers whose 
sum is < 6. Because the denominator grows quickly 
as the mode frequencies diverge, we consider corrections 


only from the two longitudinal modes closest in frequency 

to CJQQ. 



FIG. 4. (Color online) (a) Finesse as a function of length 
measured for a bare cavity. In (a-b), binned raw finesse data is 
shown through opacity, while the mean finesse as a function of 
length is indicated by the solid line, (b) Finesse as a function 
of length measured for a cavity containing a 10.5 pm diamond 
membrane. The length axis is changed from (a) so that the 
beam radius on the fiber mirror takes on the same values over 
the range of collected data, (c) A simulation of finesse as a 
function of length including perturbative coupling to higher 
order modes evaluated for different mirror clipping radii. 

To examine the importance of these perturbative cou¬ 
plings experimentally, we measured finesse as a function 
of cavity length with and without the membrane (see 
Fig. [^-b). The bare cavity finesse exhibits a decreasing 
slope as a function of cavity length, which arises from 
coupling to lossy higher order modes caused by the Gaus¬ 
sian shape of the fiber dimple [4Ql [4T] . We offset the X- 
axes of the bare cavity and membrane-in-cavity so that 
the beam radii on the fiber mirror would match, varying 
from 2.5 pm to 3.6 pm over the length range presented in 
Fig.|^-c. In addition to the decreasing slope seen for the 
bare cavity, we measured intermediate drops in finesse at 
specific cavity lengths. Since we use the same fiber mir¬ 
ror in both data sets, these dips must be associated with 
the membrane, and could be caused by coupling to lossy 
higher-order transverse modes. For comparison, we simu¬ 
lated the cavity finesse (ignoring all other loss processes) 
using the first order correction to the electric field wave- 
function (Eq. 1^ for different clipping radii on the fiber 
mirror, outside of which all light is assumed to scatter 
out of the cavity. The result is shown in Eig. |^, which 
exhibits qualitatively similar drops in finesse at certain 















resonant lengths. 

We lack the detailed surface profile data to accu¬ 
rately parameterize our membrane and fiber dimple 
topography, so the simulations cannot include the 
exact perturbations present in the measurement (for 
example, it likely also contains a wedge, which would 
couple TEMoo and TEMio modes). Nevertheless, our 
calculations demonstrate that the finesse reductions 
observed at specific cavity lengths in Eig. could 
reasonably be caused by this mechanism. Perhaps more 
importantly, this data illustrates that perturbative losses 
are not a major impediment to working with planar 
membranes, even at relatively large thicknesses > 10 
pm and over a range of cavity lengths. Eurthermore, 
as the diamond thickness diminishes, the perturbative 
coupling drops, indicating that it should be a negligible 
effect for few-micron-thick membranes. 


VI. CONCLUSION 

We have shown that high finesse ^ 17,000 can be 
maintained in a Eabry-Perot microcavity even with in¬ 
corporation of a thick diamond membrane. The mem¬ 
brane modifies the cavity modes, leading to variations 
in linewidth for different membrane thicknesses or differ¬ 
ent resonant frequencies. Our simulations indicate that 
surface losses dominate, producing qualitatively different 
behavior from bulk absorption. Membrane-induced cou¬ 
pling to higher-order transverse modes, conversely, does 
not greatly impact device performance. We anticipate 
that, despite the large surface losses, the current cavity 
will allow in the range of a 20-fold Purcell enhancement 
for diamond-based emitters, which in the case of the NV 
center would direct more than a third of its emission into 
the zero phonon line. Eurthermore, device performance 
in this case is limited by surface roughness or contam¬ 
ination that is well above the currently attainable lim¬ 
its [271[28]. 

When cooled to cryogenic temperatures and locked to 
the NV resonance frequency, such devices could signifi¬ 
cantly enhance the efficiency of photon-mediated entan¬ 
glement between distant defects [42]. Moreover, the cav¬ 
ity linewidth is below the typical spacing between spin- 
resolved resonant optical transitions in the NV center, en¬ 
abling exploration of spin-dependent cavity effects. With 
improved diamond surface treatment and higher reflec¬ 
tivity mirrors, finesse ^ 10 ^ should be possible [ 12 ], and 
shorter cavities with smaller radius-of-curvature mirrors 
could enhance the cavity cooperativity by another order 
of magnitude HZ]; even with current fabrication capa¬ 
bilities, Purcell enhancements in the range of 200 appear 
within reach. Ultimately, this highly-tunable open-cavity 
geometry could offer a route towards an efficient or even 
deterministic interface between single photons and solid- 
state spins. 
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Appendix A: Estimation of Bare Cavity Parameters 

The following section describes a method for extract¬ 
ing some parameters of the bare fiber cavity system 
given measurable quantities and a one dimensional model 
shown in Eig. These calculations permit us to under¬ 
stand the asymmetric lineshapes observed for both the 
bare and membrane-in-cavity systems. 


t t 



FIG. 5. (Color online) The one-dimensional model used to 
estimate bare cavity parameters. L is the cavity length; r and 
t are the mirror amplitude reflection and transmission coeffi¬ 
cients, Qrt is the round-trip gain of the cavity, and Ein, 

E'ef, Etrans, and Edrc arc the electric fields at the input, re¬ 
flected into the fiber, reflected into other modes, transmitted, 
and circulating in the cavity, respectively. 

In this model, the mirrors have real amplitude trans¬ 
mission and reflection coefficients t and r, and we account 
for a loss per round trip in the cavity of 1 — ^ 2a. 

Assuming that light is launched into the cavity through 
the fiber, Ein is the incident electric field, is the re¬ 
flected field that is coupled back into the fiber core, E'^^ 
is the reflected field that is not coupled into the fiber core, 
and Etrans is the field coupled to the (free space) trans¬ 
mitted mode. The field circulating within the cavity is 
represented by Edrc^ which is defined just to the right of 
the left-hand mirror; the change in amplitude and phase 
incurred in one roundtrip is represented by the rountrip 
gain: g^t = We consider imperfect cavity 

coupling, where ei is the overlap between the fiber and 
cavity modes, and €2 is the coupling coefficient between 
the cavity and transmitted modes (for our analysis we set 
62 ~ 1). If the ablation spot is not perfectly centered on 
the fiber core, the promptly reflected light that is coupled 
back into the fiber mode can be described by a complex 
coupling coefficient 77 , which has a magnitude less than 
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unity as well as a nonzero phase for imperfect alignment. 
The relevant relationships between these parameters are 
given by: 


=-e rteiEcirc + rrjEin (Al) 

-Scire = e~'^°‘~^r^Ecirc + teiEin (A 2 ) 

-^trans — ^ ^^ 2 -^circ (A3) 

Solving for the transmitted and reflected powers 
{\Etrans\‘^ l^re/P)’ normalized to the input power 
yields the power reflected (P^) and transmitted 

(Pt): 


Pr = 
Pt = 


(t2 


g4a ^4 _ 2 e2a ^2 
^2a +A ^2^2 


P44 


Aa . 


. r 2 Lu; 


C ] 


(A4) 

(A5) 


We set ?7 = a + and expand the cosine terms to second 
order in AL, where L = mA /2 + AL and m is an integer, 
resulting in power lineshapes of the form: 


Pr = 


(ai + a2 AL) (^) 


(5L\2 


(f)VAL 2 


■yo 


ns 




+ AL 2 


(A 6 ) 

(A7) 


where 5L is the FWHM cavity linewidth measured in 
length, and: 


yo = [o? + + at^e\ 


al = 'Kt^e\ 


i{r^ — + r‘^t^e\ 

|'g 2 a _ ^ 2 )^ 


a2 = 




c[e 


, 2 ck: 




a3 = TT- 


t^e\e2 


[e 


2a 


(AS) 

(A9) 

(AlO) 

(All) 


Note that this produces a Fano lineshape in reflection. 
The parameters ni, a 2 , ns} can be extracted from our 
data by fitting the transmission and reflection curves and 
calibrating the input power. In addition, we use measure¬ 
ments of the finesse P and the following relationships to 
fully constrain the cavity parameters: 

T T Of = 1, (Al3) 


as well as the known laser frequency uo and 62 = 1. With 
these expressions, one can solve for all of the cavity pa¬ 
rameters of interest. For example, using data acquired in 
a bare cavity of length 12.2 ± 0.3 /im we obtain: 


t = ( 8.8 ± 0 . 2 ) * 10 “^ 
r = -0.999957 ± 0.000001 
a = 0.61 ± 0.02 
6 = 0.14 ±0.04 
ei = 0.69 ±0.03 
a = 8 ± 1 ppm 

This yields a power transmittance of T = 78 ± 3 ppm, 
which agrees with the quoted coating value of T = 70 ± 
10 ppm. The combined absorption and scattering losses 
were quoted to be < 24 ppm, which also agrees with the 
derived a value. 

Note: While in the final stages of preparing this 
manuscript, we became aware of detailed theoretical and 
experimental exploration of the origin of asymmetric line- 
shapes associated with misalignment of fiber cavities [43] . 


Appendix B: Calculating the Effective Radius of 
Curvature 

The Gaussian-shaped ablation dimple can be approx¬ 
imated by a parabola near the center, which has a well 
defined radius of curvature. This radius is appropriate 
for cavity modes with small beam diameters on the mir¬ 
ror. As the cavity length is increased, the mode diam¬ 
eter grows and the approximation breaks down. In this 
regime, it is more accurate to estimate the effective mir¬ 
ror radius from the spacing of the higher order TEM 
modes. If Az/^rans is the difference in frequency between 
adjacent transverse modes with the same longitudinal 
mode (e.g. TEM^^^ and TEM^^^+i), the effective ra¬ 
dius of curvature is 

where L is the length of the cavity and i'fsr = c/2L is 
the free spectral range. Using this equation to analyze 
the TEMoo and TEMqi/TEMio modes in the white-light 
spectrum measured for the bare cavity, we estimate an 
effective radius of curvature of 61.0 ±1.4 qm for a bare 
cavity length of 13.3 qm and a beam radius of 2.6 qm on 
the fiber mirror. The same beam diameter would corre¬ 
spond to a cavity length of 22 qm for a cavity containing 
a 10.5 qm diamond membrane, as the beam diverges less 
in the higher refractive index medium. Deviation from 
the radius of curvature extracted from a parabolic fit to 
our interferometry measurement (L ^ 50 ± 1 qm) arises 
because the finite diameter mode samples a range of cur¬ 
vatures within the Gaussian dimple. 


Appendix C: Numerical Cavity Model 

We have developed a numerical model to calculate the 
fundamental Gaussian cavity mode for a half symmetric 
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cavity containing a diamond membrane bonded to the 
flat mirror. The model first solves for the Gaussian beam 
parameters (waist radius and position) in both the air 
and diamond sections assuming a curved diamond surface 
lying along the mode wavefront (see Fig.[^. As boundary 
conditions, we require that the beam diameters and radii 
of curvature are equal at the air-diamond interface to 
ensure electric field continuity. In addition, the radius 
of curvature in air should match the ablation radius of 
curvature at the fiber mirror, while the mode waist in 
diamond should lie at the flat mirror. With these four 
requirements, once can solve for the required radius of 
curvature for the air-diamond interface, the beam waists 
wi and W 2 corresponding to the modes in diamond and 
air, and the effective waist position X 02 for the air mode. 


Perturbation 



FIG. 6. (Color online) A diagram illustrating the 1/e^ inten¬ 
sity radius of the zero-order Gaussian modes in the diamond 
and air (solid lines). The mode in diamond has a waist with 
radius wi at the flat mirror, while the mode in air has a waist 
with radius W 2 a distance X 02 from the mirror flat. The per¬ 
turbation volume considered in Eq. (dashed lines) is the dif¬ 
ference between the presumed diamond interface lying along 
the mode wavefront and the planar diamond surface. 

Once the two Gaussian modes in the air and diamond 
regions are known, one can solve for the resonant fre¬ 
quencies and lengths of the cavity using transfer matrix 
theory applied to the left and right traveling electric fields 
within the cavity structure (see Fig. [^. 

Here, we model the full geometry of the flat mirror 
as a one-dimensional 29 layer dielectric stack where the 
right-moving electric field travels from glass into air; the 
transfer matrix describing this process is Mg^a- Gon- 
versely, the fiber mirror is modeled as the inverse ma¬ 
trix {Mamg = ^^rna)' Since the mirrors are defined for 
air termination, we need to include an additional ma¬ 
trix to model the air-diamond interface at the flat mirror 
{Dad)- As mentioned above, we assume that the sec¬ 
ond air-diamond interface {Dda) follows the curvature of 
the wavefronts at that position (dashed line in Fig. [^. 


Mirror Mirror Diamond 
Flat 


Air Mirror Fiber 


"ref 


"trans 


M D L D L M 

'•*gma ^ad "“d ^da ""a amg 


FIG. 7. (Color online) A diagram indicating the relevant 
transfer matrices used to calculate the resonant cavity fre¬ 
quencies and lengths. 


The propagation matrices in the air and diamond include 
Guoy phase, and are given by: 


L diamond 


0 ^^^^td-i(f)i{td) 


r^-^(L-td)-\-i(f)2(L-xo2) Q 

= 0 ^^{L-ta)-iML-xo2) 

where 0i(x) = arctan {xX/7rw‘f) is the Guoy phase, td is 
the diamond thickness, L is the cavity length, wi and W 2 
are the 1/e^ radii at the diamond and air waists respec¬ 
tively, and X 02 is the effective beam waist position for the 
mode in air (see Fig. |^. The full transfer matrix for the 
cavity is: 


Dfrans 

0 

where: 

S = Mamg La Dda Ld Dad Mgma- (C2) 

Using the transmission curves calculated with this model, 
the linewidth in frequency and length can be determined. 
The field within the diamond and air regions can be found 
by evaluating subsets of the transfer matrices to find the 
amplitudes in the air and diamond, which are then mul¬ 
tiplied by appropriate Gaussian modes. 

When adding loss to this model, the dielectric indices 
used in the mirror stack and diamond were given small 
complex components; in addition the interface matrices 
Dda and Dad were modified according to Eqs. [4][^ 


= 5' 


Ejr 




ref 


(Cl) 


Appendix D: Nondegenerate Perturbation Theory 

The Her mite-Gaussian family of modes represent 
the eigenstates of a spherical resonator, satisfying the 
Helmholtz equation m. Introducing some small volume 
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of material with a different refractive index can break the 
cavity symmetry, leading to a new set of cavity eigen¬ 
states that can be expressed as a linear combination of 
the unperturbed cavity modes [39] . 

In our case, the zero-order modes correspond to so¬ 
lutions in the presence of a membrane whose interface 
is curved to follow a wavefront (see Fig. [^. They are 
defined by 

(V^ + Kin^(r))?/)°(r) = 0, (Dl) 

where are the zero-order modes of the system, i^i = 
contains the corresponding eigenfrequencies and 

Uq (r) is the index of refraction inside the cavity assuming 
an air-diamond interface following the mode wavefront. 
These zero-order eigenfunctions may be found exactly as 
a (real-valued) Hermite-Gaussian family of modes with 
different parameters in the air and diamond regions (see 
Appendix The orthonormalization condition is 

jJJ (D2) 

where Smn is the Kronecker delta, the integral is taken 
over the cavity volume, and subscripts m and n encode 
all transverse and longitudinal mode indices. 

We wish to calculate the perturbative effect of the 
membrane planarity, which is equivalent to introducing 
a small piece of dieletric representing the difference be¬ 
tween the curved surface and a flat one. We are inter¬ 
ested in a fundamental transverse mode (which is non¬ 


degenerate), for which the exact eigenstate ip satisfies 

(V^ + n (^o(^) + ^^(^^))) '0 = 0^ (D3) 

where V(r) = — 1 inside the perturbation volume (and 

zero outside), A is some small number, and n corresponds 
to the new eigenfrequency. We can express both pj and 
K 8iS 8i power series in A: 


CX3 


+ 

n=l 

(D4) 

oo 

Kq + An 

n=l 

(D5) 


where and /^o correspond to the non-degenerate funda¬ 
mental mode eigenstate of the unperturbed system, and 
'0’^, An are the order corrections. Considering only 
the terms of Eq. |D3| to order A, one finds 

(V^ + nl{r)Ko)'tp^ = -(Aino(r) + KoV{r))'tpo, (D6) 

and one can thereby derive the first order correction to 
the eigenstate: 

= «„ V fffA(r)V(r)^S(r)<lr^ 

where is the mth order mode of the unperturbed 
system, and we have set A = 1. Here m labels all lon¬ 
gitudinal and transverse modes to which the zero-order 
Gaussian mode can be coupled by the perturbation. 
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